Multiple scattering of matter waves: 
an analytic model of the refractive index for atomic and molecular gases 
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We present an analytic model of the refractive index for matter waves propagating through atomic 
or molecular gases. The model, which combines a WKB treatment of the long range attraction 
with the Fraunhofer model treatment of the short range repulsion, furnishes a refractive index in 
compelling agreement with recent experiments of Jacquey et al. [Phys. Rev. Lett. 98, 240405 
(2007)] on Li atom matter waves passing through dilute noble gases. We show that the diffractive 
contribution, which arises from scattering by a two dimensional "hard core" of the potential, is 
essential for obtaining a correct imaginary part of the refractive index. 
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I. INTRODUCTION 

Multiple scattering of matter waves arises in a broad 
variety of contexts, ranging from gas transport [l| 
to the scattering of subatomic particles by atoms or 
molecules 0, @| to cosmic ray showers [4] . Inherent to the 
concept of multiple matter wave scattering is the notion 
of a refractive index, n, which, by the optical theorem 
is proportional to the forward scattering amplitude, /(O). 
While Im[/(0)] and thus Im[n] are related to the total 
integral scattering cross section and so can be measured 
in a standard beam scattering experiment, Re[/(0)] and 
thus Re[n] reflect the change of the wave's phase veloc- 
ity and can only be accessed by matter wave interfer- 
ometry. The first such interferometric measurement was 
carried out in 1995 by Schmiedmayer et al. 0,0] > who 
determined both the attenuation and the phase veloc- 
ity change of sodium atom waves propagating through 
a number of gases. Subsequently, the Pritchard group 
was able to observe glory-type oscillations in the depen- 
dence of the refractive index on the sodium beam veloc- 
ity [sl-fToj. In 2007, Jacquey et al. [fl| implemented an 
improved, Mach-Zehner type of an atom interferometer 
and measured the index of refraction of lithium waves 
passing through noble gases. In 2008, Champenois et 
al. [lH and, independently, Hornberger and Vaccini [l3| . 
provided an analysis of how the motion of the scatterers 
affects the measured refractive index. This work showed 
that the formulae used to extract the refractive index 
from the experimental data 0, and in the related the- 
oretical treatments [isl - fTsI lacked consistency with the 
Beer-Lambert law and thus were incorrect. 

According to the scenario of Ref. an atomic or 

ionic beam of "projectiles" p of mass nip and the labo- 
ratory velocity Vp propagates through a gas of density 



Nt made out of "target" particles t of mass mt and lab- 
oratory velocity v^. For plane matter waves, the veloc- 
ities Vp.t are related to the laboratory wave vectors kp^t 
by Kkp^t — 'm-p,t'^p,t- The relative motion of the target 
and projectile particles is described by the relative wave 
vector kr = tiVr/h, where = Vp — Vt is the relative 
velocity, and = mpmt/{mp + rrit) the reduced mass. 
Taking this scenario into account, Champenois et al. [l^ 
derived a new formula for the refractive index of a matter 
wave propagating through a dilute gas 
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where /(fc^, 0) = /(O) is the forward scattering amplitude 
pertaining to the scattering angle •& = 0. The averaging 
in Eq. ([T]) is carried out over the relative wave vectors, 
kr, corresponding to a normalized relative velocity dis- 
tribution [121, see below. 

Previous theoretical treatments of the refractive index 
of atomic waves, surveyed in Ref. [l^ . were based on 
WKB or eikonal [ifl approximations, which be- 
come analytic only if the repulsive, short-range interac- 
tion, is neglected, cf. also refs. [3, ilSj- In this contribu- 
tion, we combine the WKB treatment of scattering for 
the long-range attraction with the Fraunhofer approxi- 
mation for the short-range repulsion and thereby obtain 
an analytic model of the refractive index which accounts 
for both short- and long-range multiple matter- wave scat- 
tering. The resulting analytic refractive index provides 
an additional insight into multiple scattering and facili- 
tates data analysis. 

After briefly introducing the model, we apply it to the 
refractive index of Li atom waves propagating through 
Ar, Kr, and Xe gases and compare the results with the 
measurements of Jacquey et al. [ll|. In addition, we 
exemplify the model's scope by examining the refraction 
of a Na"*" ionic beam passing through a dilute N2 gas. 
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II. MODEL OF MATTER- WAVE REFRACTION then becomes, on substituting from Eq. 



In the model introduced herein, we consider the scat- 
tering amphtude to consist of two parts: 



f{kr,0) = /short(fcr,0) + /long(fcr,0) 



(2) 



The short range part, /short (fcr, 0), arises from the scat- 
tering by the "hard core," repulsive branch of the po- 
tential, and can be evaluated in closed form within the 
Fraunhofer model of matter wave scattering. The Fraun- 
hofer model was introduced by Drozdov |20j and gener- 
alized by Blair [21] in the late 1950s to treat inelastic nu- 
clear scattering, and adapted by Faubel to account 
for rotationally inelastic thermal collisions between he- 
lium atoms and N2 and CH4 molecules. Recently, we 
generalized the Fraunhofer model to treat atom-molecule 
collisions in electric [^J], magnetic [2^, and laser 
fields [26*1 , and also to gain insight into the stereodynam- 
ics of molecular collisions [13, [1^ . 

Within the Fraunhofer model, the hard core part of the 
potential amounts to a two-dimensional, impenetrable 
obstacle with sharp boundaries. The elastic scattering 
amplitude is given by the amplitude for the Fraunhofer 
diffraction of matter waves by such an obstacle [2^ , 
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where Ji is the Bessel function [2£] and Rq is the radius 
of the interaction which, for an atom-atom potential, is 
given by the solution of the equation V{Ro) = ?i^fc^/(2/i), 
with V{r) the interaction potential. For forward scatter- 
ing, the amplitude of Eq. ([3]) becomes 



/short (fcr, 0) 



(4) 



Hence the Fraunhofer amplitude is purely imaginary and, 
therefore, contributes only to the imaginary part of the 
refractive index. 

In order to account for the long-range part of the scat- 
tering amplitude, /iong(fcr,0), we make use of the WKB 
approximation, which is accurate for thermal collisions 
dominated by partial waves with large angular momenta 
I. We consider the general case of the inverse-power long- 
range potential. 



(5) 



with Cfi > and /3 > 3. For an atom-atom interaction, 
(3 — 6, while /? = 4 for an ion interacting with an atom or 
molecule. For the potential of Eq. ([S]), the WKB phase 
shift for the l-th partial wave is given by Q 
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The WKB forward scattering amplitude 
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III. REFRACTION OF ATOM MATTER WAVES 

In the case of an atom beam passing through an atomic 
gas, /3 = 6 and the forward scattering amplitude becomes 



/atom(fc., 0) = /short (fcr, 0) + //,1(A:„ 0) 



ikrRoikr) f ^kr\^^^ f StT ^Cg ^ 
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By substituting from Eq. ([9]) into Eq. ([T]) and making 
use of the distribution function given by Eq. (A3) of 
Ref. [H, 
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with a = y^2kBT/mt, ks Boltzmann's constant, and T 
the temperature, we obtain the following expression for 
the reduced index of refraction. 
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where the velocity averages are given by 



{Vr) — / VrP{Vr)dVr, 

Jo 



vl'^P{Vr)dVr. 



(12) 



(13) 



In general, the i?o value depends on the relative veloc- 
ity Vr, but for narrow velocity distributions, such as those 
implemented in the experiments of Jacquey et al. [ll[, 
the dependence is found to be negligible. We could thus 
simplify the resulting expression by fixing R^ivr) to a 
constant value pertaining to the mean relative velocity, 

i?0 = Ro{{Vr)). 

While Table |T] compares the ex per imental values of 
??atom obtained by Jacquey et al. [Tl| with our model 
results for Li matter waves propagating through Ar, Kr, 
and Xe gases. Figure 1 singles out the Li-|-Xe system and 
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FIG. 1. (Color online) Upper panel shows the real and imag- 
inary parts of the reduced refractive index r) = (n — 1)/Nt 
[m'^] multiplied by lO^^v^''^ for Li beam propagating through 
Xe gas, in dependence of the beam velocity Vp. Lower panel 
shows the corresponding ratio Re (r;)/Im {ri). Model re- 
sults (solid curves) are compared with experimental data of 
Ref. hj. Dashed curves show the results of the model with 
the short-range potential excluded. 



compares the experimental and theoretical dependence 
of the refractive index on the Li velocity. As an input 
for the analytic model we used potential energy curves of 
Ahokas et al. [sO] . The analytic refractive index is seen to 
be in a compelling agreement with the experiment. We 
note that the interference of the matter waves scattered 
by the short- and long-range parts of the interatomic po- 
tential causes the refractive index to exhibit glory-type 
oscillations. These are not correctly rendered by our an- 
alytic model, as the approximate scattering amplitudes, 
which were obtained by different methods, fail to "keep 
track" of one another. 

Also presented in both Table H] and Figure 1 are the 
values of the refractive index for the long range potential 
included but the short-range interaction excluded, i.e., 
for /(fcr,0) = /iong(fcr,0), cf. Eq. ([2]). One can see that 
in such a case the imaginary part of the refractive index 
is substantially less than in the experiment. 

In the limit of a cold atomic gas, i.e., for T — > 0, the 
velocity distribution (ITUl) becomes a (5-function, which 
simplifies Eq. PT|) . since then \ ^ " o^^^ 
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— >■ Vp and 
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In this case, Re[r7atom] is proportional to Vj 
while Im[77atom] is proportional to both and Vp 
However, since the term arises from the diffraction 
of matter waves from the hard core of the potential, it 
dominates Im [ryatom] at large collision energies at the 



expense of the Vp 
range attraction. 
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term, which arises from the long- 



IV. REFRACTION OF ION MATTER WAVES 



TABLE I. Reduced index of refraction, rj = [n ~ 1)/Nt [m^], 
for a Li beam propagating through different noble gases with 
a mean velocity Vp = 1075 m/s. 





Experiment 
of Ref. [11] 


Model 


Model 
/short (A:r,0) = 


Ar 


10^^Re(»7) 


1.20 ±0.11 


0.97 


0.97 


1029lm(j7) 


2.11 ±0.06 


1.65 


1.33 


Re(j7)/Im(r,) 


0.56 ±0.05 


0.58 


0.73 


Kr 


10^^Re(r7) 


1.57 ±0.10 


1.25 


1.25 


102''lm(r;) 


1.99 ±0.07 


2.04 


1.71 


Re(?7)/Im(r7) 


0.78 ± 0.04 


0.61 


0.73 


Xe 




1.82 ±0.07 


1.52 


1.52 


102«Im(?7) 


2.40 ± 0.06 


2.42 


2.08 


Re(7))/Im(77) 


0.70 ± 0.03 


0.62 


0.72 



For an ion beam passing through an atomic or molec- 
ular gas, /3 = 4 and the scattering amplitude becomes 



/ion(fc., 0) = /.hort(fcr, 0) + fi^^^Ukr, 0) 



(14) 



2 V 2 / V3, 

leading to a reduced refractive index for an ionic beam: 

n — 1 



27r 



i{vr)Rl 
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with 



,l/3\ 



(15) 



(16) 



We note that in the case of an anisotropic interaction 
of, say, an ion or atom with a linear molecule given by 
a potential V{r,9), the Rq value is given by the spheri- 
cal part of the hard core of the potential [23j. In order 
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TABLE II. Reduced index of refraction, ri={n- l)/Nt [m^], 
for a Na""" beam propagating through N2 gas with a mean 
velocity Vp = 5000 m/s. 





Model 


Model 






./short (fcr, 0) = 


1029Re(?7) 


0.25 


0.25 


102''lm(r?) 


0.16 


0.14 


Re(7;)/Im(r,) 


1.60 


1.73 



to extract Rq, we first solve the equation V{R{9),9) — 
/i^A;^/(2/i) to obtain the hard core shape R{0), which we 
then expand in a Legendre series, 

R{e) = i?oPo(cos6') + i?iPi(cos6i) + R2P2{cos0) + 

(17) 

from which we deduce the Ro value. 

The values of the refractive index for an ion beam 
passing through a molecular gas, as exemplified by the 
Na^-N2 system, are listed in Table lU at an ion projec- 
tile velocity Vp = 5000 m/s. The requisite value of R^ 
was extracted from the potential energy surface of Ref. 
[3l| . We see that the refractive index is less affected by 
switching off the repulsive part of the potential for the 
ion-molecule system than for the atom-atom system, cf. 
section IIIII This is due to the greater strength of the 
—C4r~^ interaction compared with that of the —Cer~^ 
potential and the more dominant role it thus plays in 
determining n. 

We note that in the limit of zero temperature, Vr Vp, 



the real part of the refractive index (ITSt becomes pro- 
portional to Vp . The imaginary part of the refractive 
index is determined by two terms, one of which is pro- 
portional to and the other to 



-5/3 



V. CONCLUSIONS 

We presented an analytic model of the refraction of 
atom or ion matter waves passing through atomic or 
molecular gases. The values of the refractive index fur- 
nished by the model were found to be in good agreement 
with experiments of Jacquey et al. 11| . Our analysis has 



shown that in order to appraise the imaginary part of 
the refractive index correctly, we need to account for the 
diffraction by the repulsive hard core part of the inter- 
action potential. This can be achieved by invoking the 
Fraunhofer model. 

Within the Fraunhofer model of matter wave scatter- 
ing, the elastic scattering amplitude depends solely on 
the spherical radius, Rq, of the scatterer. As a result, an 
external field (e.g., magnetic), does not affect the elastic 
cross section furnished by the model. The corollary is 
that the field would not affect the index of refraction 
obtained from the analytic model presented herein either. 
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